Fluid pumped by magnetic stress 
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A magnetic field rotating on the free surface of a ferrofiuid layer is shown to induce considerable 
fluid motion towards the direction the field is rolling. The measured flow velocity i) increases with 
the square of the magnetic field amplitude, ii) is proportional to the thickness of the fluid layer, and 
iii) has a maximum at a driving frequency of about 3 kHz. The pumping speed can be estimated 
with a two-dimensional flow model. 

PACS numbers: 47.65.+a, 75.50.Mm, 47.62. +q, 47.60.+i 



Polarizable fluids can show a macroscopic reaction to 
external electric or magnetic fields. While for most fluids 
the influence of a magnetic field is fairly weak, colloidal 
suspensions of magnetic particles - so called Ferrofluids - 
do show a strong response particularly to static magnetic 
fields 0. If these fields are time dependent, a rich variety 
of phenomena occurs. The internal rotation of the mag- 
netization in an externally rotating magnetic field gives 
rise to nontrivial effects, as summarized in a recent spe- 
cial issue |2j. A particularly interesting example is the 
driving of a macroscopic flow by means of an external 
rotating magnetic field (3, because it should allow for a 
very fine tuning both of the speed and of the direction of 
the flow even in microscopic channels p|. 

In this paper we present a novel magnetic fluid p] in an 
open flow geometry especially designed for a quantitative 
comparison between the measured flow velocity and its 
theoretical estimation. It is scetched in Fig. 1. 




FIG. 1: Experimental setup. The arrows indicate the direc- 
tion of the oscillating magnetic fields provided by the coils. 

A circular Macrolon® duct with a mean diameter d of 
100 mm and a square cross-section of 5 mm x 5 mm (and 
2.5 mm x 2.5 mm in a second set of measurements) is 
filled brimful with a magnetic fluid fU with a viscosity 
of i] = 5.4 • I0~ 3 Pas. The orientation of the two coils 
producing the rotating magnetic field is also indicated in 
Fig. 1: One coil is wrapped around this circular channel 
and provides a magnetic field in azimuthal direction, and 
the outer coil produces the vertical component of the 
field. Both coils are driven with an ac-current with a 
phase difference of 90°, thus producing a rotating field 



on the free surface of the fluid within the duct. 

The characterization of the magnetic susceptibility of 
the fluid in the duct is of primary importance for the 
pumping of the fluid, as discussed in this paper. It has 
been measured as a function of frequency of the exter- 
nal oscillating azimuthal magnetic field by means of a 
pick-up coil placed into the liquid. Precisely speaking, 
the magnetization was determined from the difference of 
the signal detected by the pick-up coil in the empty and 
the filled channel, under the influence of an oscillatory 
azimuthal field. The results are presented in Fig. 2. The 
measured susceptibilities of this novel cobalt based fluids 
are fairly large compared to the more common magnetite 
based fluids. In particular, the large imaginary part of 
the susceptibility is essential for the pumping action de- 
scribed in this paper. 
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FIG. 2: The magnetic susceptibility as a function of the fre- 
quency of the external alternating magnetic field. Here the 
real and imaginary parts are denoted by squares and circles, 
respectively. 

The pump does work: A rotating field produced by the 
coils leads to a motion of the fluid in azimuthal direction 
of the channel |fj. Its velocity is on the order of mm/s 
and can thus be determined by visual inspection of tracer 
particles swimming on its surface. The next observation 
is also a qualitative one: When changing the phase dif- 
ference between the ac-fields from +90° to —90°, the flow 
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changes its sign. The flow direction is such that the vor- 
ticity of the flow field is locally parallel to the rotation 
vector of the magnetic field, i. e. the fluid flows towards 
the direction the field is rolling. 

For quantitative measurements it is necessary to use 
particles which are small compared to the channel width 
(dandruff, diameter about 1 mm). The velocity is de- 
termined by taking the time a particle needs to travel 
a few centimeters. The size dependence of these mea- 
surements is taken into account by using the numerically 
calcultated - roughly parabolic - velocity profiles , and 
by assuming that a floating particle represents the mean 
speed with respect to its diameter. A result obtained for 
a fixed frequency of 1 kHz is presented in Fig. 3, where 
the maximal flow velocity within the channel is shown 
as a function of the amplitude Go of the driving exter- 
nal magnetic field. The velocity increases proportional to 
Gg, as demonstrated by the solid line, a parabola. Here 
we follow the nomenclature of Ref . , where the external 
magnetic far field is marked by G and the local one by 
H. 
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FIG. 3: Maximal velocity measured as a function of the field 
amplitude at a fixed frequency of 1 kHz. 

Having demonstrated this quadratic dependence of the 
velocity on the magnetic field, a first approach to collapse 
data obtained at different fields is to introduce a reduced 
velocity by dividing with the square of the external field. 
Another important influence determining the fluid veloc- 
ity is the height of the channel: It turns out that the 
velocity is larger in bigger channels. This leads us to re- 
duce the velocity also by dividing by the height L of the 
duct. In order to get a dimensionless number one also has 
to scale with the viscosity of the fluid. Thus we define 
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as a reduced flow velocity. Its measured values are pre- 
sented as a function of the driving frequency of the ro- 
tating magnetic field in Fig. 4, for velocities obtained in 
a 5 mm x 5 mm and a 2.5 mm x 2.5 mm duct. Both 
measurements show a maximum of this velocity in the 
range of 2 - 4 kHz. 
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FIG. 4: The reduced velocity u as a function of the driving 
frequency. Solid circles (open squares) are obtained in the 5 
mm x 5 mm (2.5 mm x 2.5 mm) duct, respectively. The 
solid line represents the value expected on the basis of the 
measurement of the ac-susceptibility. 



The pumping action can be understood as a manifes- 
tation of the magnetic stress acting on the magnetized 
fluid, as summarized in Ref. 8J. This assumption ex- 
plains all qualitative features of the observation: as long 
as the magnetization is proportional to the magnetic field 
(which has been measured to be the case for our fluid, 
within a precision of 5 % for fields up to about 1500 
Am -1 ), the stress must be proportional to G 2 as demon- 
strated in Fig. 3. If the frequency approaches zero, the 
magnetization and the field are parallel to each other, the 
tangential stress is zero and thus the motion of the fluid 
stops. For finite frequencies the velocity is proportional 
to the x ,f component of the susceptibility, which must in- 
crease linearly (to lowest order) with the frequency. For 
higher frequencies, the imaginary part of the susceptibil- 
ity %// has a maximum at about 1 kHz, which explains 
that the maximal pumping velocity is observed around 
that frequency. 

For the quantitative calculation we solve the Laplace 
equation for the velocity field including the no-slip con- 
dition for the fluid at the bottom and the side walls of 
the duct. At the upper surface of the fluid the magnetic 
stress provides the boundary condition 



dv no , 
VlT = ^(M Z H X -M X H Z 
oz 2 



(2) 



because Eq. (65) of Ref. |8J is applicable to our geometry. 

Both field components H x and H z and the correspond- 
ing magnetization M x and M z are not constant in our 
case, but depend both on time and space. By numeri- 
cal computations of the internal field H and the result- 
ing magnetization M — xif, which will be published 
elsewhere, we get for the maximum fluid velocity in the 
middle of the upper surface 
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Here N e g denotes an effective demagnetization factor and 
a a reduction factor due to the given geometry depending 
on the aspect ratio a = L y /L z of the rectangular channel. 
In our case a — 1 and we obtain a w 0.369. 

In our square cross section, the magnetization - and 
the ensuing stress at the surface of the fluid - is not ho- 
mogenous. We thus extract an effective demagnetization 
factor iVeff by fitting Eq. J2J| to the numerically obtained 
velocity for different values of X- We obtain N e s ~ 0.656, 
which seems realistic when comparing to N = 0.5 for the 
case of a circular cylinder. 

From Eqs. 111.'?!) we finally get the theoretical estima- 
tion for the reduced velocity 

a + jVoffX') m 

U 2(1 + N cSX iy + {N cSX ") 2 [ ' 

The limiting case of an infinitly wide channel (N = 1, a = 
1) is included in this formula. 

The reduced velocity obtained from Eq. is presented 
in Fig. 4 as a solid line, with the values of x 1 and x 11 taken 
from the measurements presented in Fig. 2. The agree- 
ment between the measured velocities and the solid line 
is on a 20% level. This can partly be attributed to the 
limited accuracy of the measurement procedure, which 
is indicated by the error bars in Fig. 4. The systematic 



deviations between the data and the theoretical curve 
are believed to reflect the precision of the simplifying as- 
sumptions going into the consideration presented above. 
For example, the influence of the shape of the meniscus 
of the fluid adds a three-dimensional complication to the 
problem, whose influence on the maximal pumping speed 
is hard to estimate. Moreover, it should be noted that 
magnetic fluids are not perfectly stable both in their mag- 
netic susceptibility and their viscosity, which might add 
to the small mismatch between the expectation based on 
the measurement of the susceptibility and the measured 
velocity. Finally, taking into account the small ampli- 
tude of the magnetic field, any magnetoviscous effects 
have been neglected for the calculations. 

In summary, the pump presented here works well and 
has an interesting potential especially in small geometries 
where a mechanical driving of the flow is not possible. 
More importantly, it does seem safe to conclude that the 
ansatz of a magnetic stress driven motion captures the 
essence of this pump on a quantitative level. 
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